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Current and future experimental studies of tlie liigli-energy limit of QCD, dominated by non- 
^ . linear gluon dynamics in the low- a:: sector of the hadron wavefunctions, are presented. Results 

' at HERA (proton) and RHIC (nucleus) pointing to the possible onset of parton saturation 

' , phenomena, and perspectives at the LHC and new proposed DIS facilities are outlined. 
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J> ' 1 Introduction 

(N 

Gluons provide the dominant contribution to hadronic scattering cross sections at high-energies. 
Deep-Inelastic (DIS) experiments of electrons off protons at the HERA collider at DESY have 
\Q . shown that for values of the parton momentum fraction x = PpaHon/p proton ^ 0.01, the proton 

O \ wavefunction is basically purely gluonic (Fig. [1]). This is so, because gluons are "cheap" to 

^ ■ radiate: the probability of emitting a gluon increases as oc ag \ii{Cp) and oc as \n{l/x) according 

' to the standard linear QCD evolution (DGLAP^^ and BFKlI2] resp.) equations. As a matter 

of fact, for decreasing values of x the gluon density increases so fast that unitarity would be 
^ \ ultimately violated, even for processes with large virtualities ^> Ag^j-,. The theoretical 

\ expect at ionlSHl is that at some small enough value of x (a<j ln(l/x) ^ 1) non-linear gluon-gluon 

fusion effects - not accounted for in the DGLAP/BFKL equations ~ will become important and 
will tame the growth of the parton densities. The onset of saturation in the proton (or in a 
nucleus with A nucleons) is expected for parton momenta < where Qs is a dynamical 
"saturation scale" ^^El -which depends on the transverse size (vri?^) of the hadron: 

Ql{x) ^ as-^ xG{x, Q2) ^ ^1/3 ^-A ^ ^i/3(^)A ^ A^/^e^\ (1) 

with A ~ 0. 250. Eq. ([1]) tell us that Qs grows with the energy of the collision, ^/s, or equivalently, 
with the rapidity of the parton y = ln(l/x). The nucleon number A dependence implies that, at 
equivalent energies, saturation effects will be amplified by factors as large as A^/^ ~ 6 in heavy 
nuclear targets {A = 208 for Pb) compare d to protons. The regime of high gluon densities is 
often described in terms of the colour-dipole^SE! or "Colour Glas s Condensate" (CGC)I^ effective 
theories, with the corresponding non-linear 

bk/jimwekEHEI 

evolution equations. 



2 Gluon saturation at HERA ? 



Although the arguments for saturation are well justified theoretically, no strong deviation from 
the linear QCD equations has been conclusively observed in the perturbative kinematical range 




Figure 1: Left: Structure function F2{x,Q^) measured in proton DIS at HERA (y^s — 320 GeV) and fixed-target 
(y^ « 10-30 GeV) experiments. Right: Gluon distribution function in the proton, xG{x,Q^ = 5 GeV'^), derived 
from DGLAP-based analyses of F2 scaling violations (for comparison, the bottom curve shows the u quark PDF). 



covered at HERA. Most of the experimental observables, in particular those of more inclusive 
nature such as the total a-y*p cross-section d'^a/dxdQ'^ Itto^ / {x Q^)F2{x,Q'^), are in good 
accord with the standard DGLAP expectations (Fig. [1]). However, it is worth to note that the 
saturation scale at HERA energies is in a regime of relatively low virtualities, ~ 1 GeV^ 
and thus - since Bjorken x and virtuality are correlated as a; ~ Q"^ / s - the interpretation of 
most of the truly low-x range probed (x < 10"'^) is "blurred" by its proximity to the non- 
perturbative regime. As a consequence, the gluon distribution function xg{x,Q'^) indirectly 
obtained from the F2 scaling violations - via dF2{x,Q'^)/dln{Q'^) ^ 10 as{Q'^)/{27TT) xg{x,Q'^) 
- is poorly constrained below x ~ 10~^. Different DGLAP parametrization^l based on DIS-only 
data (ZEUS-PDF, Hl-PDF, Alekhin02) or on global-fits (CTEQ6.5M, MRST-NLO) yield gluon 
PDFs differing by factors of 3 or more (Fig. [U right). 

Notwithstanding those uncertainties, there are three empirical observations at HERA (summa- 
rized in Fig. [2]) that favour a possible onset of low-x parton saturation in the proton. The 
leftmost plot shows the geometric scaling ' property of inclusive o"dis which, instead of being 
a function of x and separately, for x < 0.01 it features a single dependence on the param- 
eter r = Q'^/Qlix) where Qs{x) = Qq{x/xq)^ with A ~ 0.3, Qq = I GeV, and xq ~ 3-10-'^. 
Such a scaling property is naturally explained by gluon saturation models, whereas the DGLAP 
approach can only reproduce it via a fine tuning of the initial parameterization of the gluon 
distribution used. Another piece of evidence for saturation effects at HERA is provided by 
diffractive processes, where the proton remains intact after the "quasi-elastic" interaction with 
the photon. Diffractive scattering, accounting for 10-15% of the total DIS cross-section, is char- 
acterized by colourless two-gluon exchange, and thus it constitutes a sensitive probe of the gluon 
densities^l^. Surprisingly, the ratio of the diffractive to total cT^*p cross-sections is found^l^ to 
be roughly constant as a function of the center-of-mass energy W and (Fig. [2l center). 
This is in disagreement with the naive pQCD expectations of an increase of the ratio according 
to rf^f = adiff/o-tot oc \xg{x,Q'^)\'^ /xglx^Q"^) ~ W^^ /W"^^ ~ W'^^. The last indication of a 



"See http : / / durpdg . dur . ac . uk/hepdata/pdf 3 . html 



"tension" between the standard linear QCD equations and low-x HERA data comes from the 
longitudinal structure function, Fi{x, Q"^) which at variance with i<2, is directly proportional to 



xg{x, The Fl{x, Q^) derived from NLO DGLAP analyses becomes unphysically negative^ 
below X ~ 10~^ for relatively small values, whereas it is a well-behaved object in saturation 
models (Fig. [2l right). 
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Figure 2: Hints of saturation at HERA. Left: Geometric scaling''' of as a function of r = QVQ? for x < 0.01. 

12 2 
Center: Ratio of diffractive to total DIS cross-sections as a function of W and Q . Right: L ongitudinal structure 



function, Fl{x,Q^= 2 GeV^), predicted by DGLAP^^ and saturation' 
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models. 



3 Gluon saturation at RHIC ? 

Among the interesting observations in nucleus-nucleus (A-A) collisions at RHIC is the possible 
onset of parton saturation phenomena. Though nuclei at RHIC are probed at lower energies 
i\/^NN ~ GeV) than protons at HERA, saturation effects are "amplified" thanks to the 
increased transverse parton density in the former compared to the latter. Two empirical obser- 
vations support the Color-Glass-Condensate (CGC) predictions of a reduced parton flux in the 
incoming ion s due to enh anced non-linear QCD effects. On the one hand, the measured hadron 
multiplicities flSIIZIIHIISl^ dNch/drjL=o « 700, are significantly lower than the dN ch / dri\r^=Q ^ 

1000 values predicted by minijet^^ or Regge^^ models, but are well reproduced by CGC ap- 
proaches'^Sl. Parton multiplicity distributions at high energies are perturbatively calculable in 
saturation approaches, since they are governed by a semihard saturation scale oc y^s^^ 

with an exponent A constrained by e-A data . Assuming parton-hadron duality, hadron multi- 
plicities at mid-rapidity rise proportionally to times the transverse (overlap) area, a feature 
that accounts naturally for the experimentally observed factorization of ^ys^J^- and centrality- 
dependences in dN^h / dr]\f^=Q (Fig. O left). 

The second manifestation of CGC-like eff'ects in the RHIC data is the BRAHMS observational^ 
of suppressed yields of semi-hard hadrons {pT ~ 2 — 4 GeV/c) in d-Au relative to p-p collisions 
at increasingly forward rapidities (up to « 3.2, Fig. [3l right). Hadron production at such small 
angles is theoretically sensitive to partons in the Au nucleus with x™" = [pt/ ^/s^) exp(— r/) « 
C(10~^)^. The observed nuclear modification factor, RdAu ~ 0.8, cannot be reproduced by 
pQCD calculationsl23IlSl26]that include the same leading-twist nuclear shadowing that describes 
the d-Au data at y = 0, but can be des cribed by CGC approaches that parametrise the Au 
nucleus saturated gluon wavefunction^SZES^ 




Figure 3: Hints of saturation at RHIC. Left: Normalized dNch / dri\rj=o as a function of cm. energy and centrality 
(given in terms of the number of nucleons participating in the collision, A'part) measured by PHOBOS in Au-Au' — ' 
compared with saturation predictions'—^. Right: Nuclear modification factor RdAu ipr) for n egativ e hadrons at 
■n = 3.2 in d-Au at w's];;^ = 200 GeV: BRAHMS data^^ compared to pQCdI^^I^ and predictions. 



4 Low-x QCD studies at the LHC 



The Large Hadron Colhder (LHC) at CERN will provide p-p, p-A and A- A collisions at ^/s^ = 
14, 8.8 and 5.5 TeV with luminosities £ ~ 10^'*, 10^^ and 5-10^^ cm~^ respectively. Following 
Eq. ([l]), the relevance of low-j; QCD effects will be significantly enhanced due to the inc reased: 
center-of-mass energy, nuclear radius (A^^^), and rapidity of the produced partonsl^SESI. At the 
LHC, the saturation momentum ~ 1 GeV^ (proton) - 5 GeV^ (Pb) will be more clearly in 
the perturbative regime'^, hard probes will be copiously produced, and the x values experimen- 
tally accessible will be much lower than at previous colliders: X2 ~ 10~^(10~^) at central (very 
forward) rapidities. All LHC experiments have interesting detection capabilities in t he forw ard 
direction which will help to constrain the PDFs in the very low-x regime: (i) can 
measure inclusive jet and Drell-Yan production down to x ~ 10~^ using the CASTOR calorime- 
ter at rapi dities 5.5 < rj < 6.6 as well as Mueller-Navelet dijets (very sensitive to non-DGLAP 
evolution) S3]34| separated by rapidities as large as Ar] ~ 10, (ii) ALICE and LHCb feature 
a forward muon spectrometer (covering 2<- rj < 5) which gives them access to heavy-quark, 
quarkonia and gauge boson measurements'^ down to x ~ 10~^. 

The advance in the study of low-x QCD phenomena will be specially substantial in Pb-Pb colli- 
sions. In saturation models, there is a one-to-one correspondence between the effects of rapidity- 
and ^s^j^ -dependences, because a parton distribution boosted to higher rapidity y is equivalent 
to a distribution sampled in a process at higher ^/s^. As a consequence, the saturation physics 
explanation of the increasing forward suppression of semi-hard yields at RHIC will imply a sig- 
nificant A-A hadron suppression at LHC mid-rapidities. CGC predictions for charged hadron 
multiplicities in central Pb-Pb at 5.5 are dN / dr]\r^=Q ^ 1500, i.e. 3-4 times lower than 

the pre-RHIC era results. [As a matter of fact, such low multiplicities help CMS'^ and AT- 
L As'^ to become very competitive experiments in the heavy- ion running mode] . Arguably, one 
of the cleanest way to study the low-x structure of the Pb nucleus at the LHC is via Ultraperiph- 
eral collisions (UPCs)'29| in which the strong electromagnetic fields (equivalent flux of quasi-real 
photons) generated by the colliding nuclei can be used for photoproduction studies at maxi- 
mum energies ^/s~^ ~ 1 TeV, 3-4 times larger than at HERA. Full simulation-l-reconstruction 



studiesl2332] of quarkonia photoproduction (7Pb^ T Pb) tagged with very-forward neutrons, 
show that CMS can carry out detailed pt,V measurements in the dielectron and dimuon T decay 
channels. Such processes probe x values in the nucleus as low as x ~ 10~^ (Fig- HI left). 




Figure 4: Kinematic {x,Q^) plane probed in e-,7-A processes: existing data compared to (a) UPC T photopro- 
duction processes (left), and (b) new proposed nuclear DIS facilities (right): LHeC and EIC/eRHIC. 



5 Proposed future deep-inelastic facilities 

Two different collider projects - the Large Hadron Electron Collider (LHeC)'^and the Electron 
Ion Collider (EIC)I2I- have been recently proposed to study deep-inelastic lepton-hadron (e-p, 
e-d and e-A) scattering for momentum transfers as large as 10^ GeV^ and for Bjorken 
X down to the 10^^. Both projects have a strong focus on the study of the low-x gluon 
structure of protons and nuclei and on non- linear QCD evolution. LHeC (EIC/eRHIC) pro- 
poses to add an extra = 70-GeV (10-GeV) electron ring to the LHC (RHIC/JLab) pro- 
ton/nucleus collider(s). In the nuclear DIS sector, LHeC (EIC/eRHIC) would allow e-A colli- 
sions at ^/s^ = 2-sjE(. = 880 (63) GeV. Fig. fright, shows the {x,Q'^) ranges accessible 
to both machines. Both proposed facilities would significantly extend the (meager) kinematical 
regime of the existing nuclear DIS data, fully mapping out the range of Bjorken-x at virtualities 
around the saturation scale (Q^ 1 - 10 GeV^, indicated by the black curve in the plot) and 
providing very valuable insights on the high-energy limit of QCD. 
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